Transient tunneling current of single electron transistors 
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The transient tunneling current of single electron transis- 
tors (SETs) is theoretically investigated. The time-dependent 
current formula given by Jauho, Wingreen and Meir [Phys. 
Rev. B 50, 5528 (1994)] is applied to study the temperature 
effect on the transient current through a single quantum dot 
embedded into asymmetry barrier. It is found that the tunnel- 
ing rate ratio significantly influences the feature of transient 
current. Finally, the oscillation structures on the exponential 
growth transient current of single hole transistors composed of 
germanium quantum dots is analyzed. 



I. INTRODUCTION 

The transport properties of single electron transistors 
(SETs) have been extensively studied theoretically and 
experimentally. 1-6 The main structure of SETs consists of 
a single quantum dot (QD) and three electrodes (source, 
drain and gate). The manipulation of SETs is based on 
Coulomb blockade effect arising from the particle interac- 
tions of QD. From the practical point of view, it is impor- 
tant for SETs to operate at room temperature. Therefore, 
the size of QD is required to be less than 10 nm for silicon 
(Si) or germanium (Ge) semiconductor QDs. In such size 
range of QD, the quantum confinement effects and charg- 
ing energies of QDs are larger than the thermal energy of 
room temperature. Recently, the Coulomb oscillation and 
staircase features of tunneling current of room tempera- 
ture SETs have been reported by several groups, where Si 
or Ge QDs are embedded into S1O2 matrix 4-6 . 

Even though it is difficult to align a single QD of 
nanometer with electrodes in the fabrication of individual 
$ETs, several methods were used to solve this difficulty. 
Nevertheless, those methods still can not precisely con- 
trol the barrier thickness, which significantly influences 
the tunneling time for electrons to access QDs. Further 
understanding for the location of QDs, the implementa- 
tion technology of SETs using thermal oxidation method 6 
can be improved. The measurement of transient current 
can provide above information. In this study we apply 
the current formula of Ref. [2] to examine the tran- 
sient current of the SET with asymmetrical tunneling 
rates at room temperature. It is found that the tran- 
sient current exhibits exponential decay when electrons 
are injected into the QDs from the left electrode and 
Tl/Tr < 1, where T L / R represents the tunneling rate 
for electrons from the left (right) electrode to the QD, on 



the other hand the transient current exhibits the expo- 
nential growth as Tl/Tr > 1. In addition to theoretical 
studies of transient current, we report the experimental 
measurement of transient current of single hole transistor 
(SHT). 



II. FORMALISM 

The transient current of a SET consisted of a sin- 
gle QD and three electrodes was theoretically derived 
by Jauho, Wingreen and Meir. 2 When applied voltage 
is insufficient to overcome the charging energies aris- 
ing from electron-electron repulsion interactions in the 
QD,the time-dependent tunneling current through the 
ground state of QD can be expressed as (ref. [2]) 
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According to Eq. (1), there are four components for the 
net current from the left electrode to the right electrode. 
The current of Eq. (2), J%j R (t), denotes the current flow- 
ing out from the QD to the left (right) electrode. This 
current results from the time-dependent electron occupa- 
tion number of QD, which is given by 



N(t) 



2tt 



(T L f L (e)+T R f R (e))\A c (e,t)\ 



(4) 



J L n / R (i) denotes the current flowing into the QD from the 
left (right) electrode. Obviously, J™/ R it) is determined by 
the time-dependent spectrum function ImAc(e,t) where 
Im means taking the imaginary part of Ac{e,t). Nota- 
tions e and Ti denote, respectively, the electron charge 
and Plank's constant. For the simplicity, we assume 
that the tunneling rates T^/ R are bias and energy in- 
dependent. Ft, and T R denote, respectively, the tunnel- 
ing rates from the left and right electrodes to the QD. 
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fL(R)(e) = l/(exp^-' lL w'>/ kBT + 1) is the Fermi distri- 
bution function of the left (right) electrode. The chemi- 
cal potential difference is related by the applied voltage 
Ml — = eV a - From the results of Eqs. (3) and (4), the 
time-dependent spectrum function plays a crucial role in 
the determination of time-dependent tunneling current. 
The system studied here considers that the step- like mod- 
ulation is added into the gate electrode, the expression of 
Ac(e,t) is obtained as 



A c {e,t) = 



1 



e- (E 1 -eV g )+iY/2 
1 



(5) 



L e- (E 1 - eV g ) +iT/2 e-E^iT/2' 

where E\ is the ground state energy level of QD. T = Tl + 
Tn. In the absence of gate voltage V g the Ac(e) becomes 
time- independent retarded Green's function A c (e) = 
1/ (e — E\ + iY /2). The resonant energy level Ei is shifted 
to E 1 -eV g , that is A c (e) = A c (e) = l/(e- (E 1 - eV g ) + 
iT/2), when the system goes into steady state. In the 
transient process the density of states of QD depends on 
time. Consequently, J(t) displays time-dependent behav- 
ior in the transient process. 



III. RESULTS AND DISCUSSION 

Due to the complicate spectrum function, the transient 
current lacks analytic form. To numerically calculate J(t) , 
we set the Fermi energy level of electrodes and the ground 
state energy level Ep = 50 meV and E\ — 90 meV, re- 
spectively. Therefore, the energy levels of QD are empty 
at zero temperature under zero bias. First of all, we plot 
the time-dependent tunneling current through a symmet- 
ric double-barrier tunneling structure with Tp = 0.5 meV 
and = 0.5 meV in response to a step-like modulation 
for different applied voltages at zero temperature in Fig. 
1: solid line (V a = 30 mV), dashed line (V a = 25 mV) 
and dotted line (V a = 20 mV). The tunneling current is 
zero at t — since the Fermi level of electrodes is below 
the resonant state E\. When the gate voltage is added 
into the system, electrons from the left electrode are in- 
jected into the new resonant energy level E\ — eV g . Conse- 
quently, the tunneling current jumps instantly and finally 
reaches the steady state. Current also displays an inter- 
esting ringing behavior, which was pointed out in Ref.[2]. 
The oscillations of the curve do not maintain constant 
frequency in the transient process. In addition, the fre- 
quency of oscillations is increased, when the separation 
between the Fermi energy of left electrode and the reso- 
nant level Ei is decreased. Due to symmetry tunneling 



rates, the behavior of J(t) can be understood by the anal- 
ysis of J™(t) = J£M) + J™Jt) - J™ (t), where 
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We plot J^t) (dashed line) and Jfy(t) (dotted line) in 
Fig. 2. Jp\(t) exhibits the behavior of (1 — exp~^\ but 
J L ™2 W exhibits the oscillation. Owing to the Fermi energy 
of right electrode away from E\ — eV g , we see small con- 
tribution from the J^(t). From the results of Fig. 2, the 
oscillation of the solid curve shown in Fig. 2 is consisted of 
two frequencies u)\ and uo 2l which are, respectively, given 
by J™ 2 (t) and J™{t). 

Eq. (7) reveals that the frequency of oscillation de- 
pends on not only V a , but also V g . To demonstrate 
it, we show J(t) for different applied gate voltages and 
V a = 30 mV in Fig. 3: solid line (V g = 25 mV), dashed 
line (V g = 20 mV), dotted line (V g = 15 mV) and dot- 
dashed line (V g = 10 mV). We see that the frequency of 
oscillation can be tuned by the magnitude of applied gate 
voltage. For V g = 10 mV, the Fermi energy level is just 
aligned with E\ — eV g , the oscillatory behavior of tunnel- 
ing current almost vanishes. Although the charge density 
of QD is not be directly measured, we plot N(i) in Fig. 
4, where the curves correspond to those of Fig. 3. It is 
very clear that QD is empty before the gate voltage turns 
on. When the system reaches the steady state, the value 
of N(t) is fractional due to the result of opened system. 
It is unexpect that N(t) does not exhibit very manifest 
oscillation feature as well as J(i). 

Because the width of tunneling barrier could not be pre- 
cisely controlled in the fabrication of SETs, we attempt to 
study the tunneling current for the case of asymmetrical 
tunneling rates. We show J(t) for different tunneling rate 
ratio in Fig. 5: the curves from the bottom to the top 
(for t/t < 0.5) correspond to T L = 0.3,0.4,0.5,0.6 and 
0.7 meV, respectively. Meanwhile T = Tl+Tb = 1 meV. 
Note that other parameters used in Fig. 5 are the same as 
those of the solid line of Fig. 1 . It is worthy noting that 
the maximum steady current occurs for the symmetry 
tunneling rates. This is consistent with other theoretical 
reports 1 . In addition, the steady current of Tl = 0.7 meV 
(Tl = 0.6 meV) is the same as that of Tl = 0.3 meV 
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(r^ = 0.4 meV). This indicates that the formula given 
by Eq. (1) maintains the asymmetry of Tl and Tr. 

It is not avoided for high temperature SETs to study 
the temperature effect on the transient current. Figure 
6 shows J(t) with symmetry tunneling rates for different 
temperatures at V a = 30 mV and V g — 25 mV. The 
ringing behavior of J(t) is swept out by temperatures. 
In addition, the magnitude of J(t) is suppressed with in- 
creasing temperature since the right electrode provides 
current J R ut (t) into the QD, which is opposite to J™(t). 
Consequently, the net current becomes small with increas- 
ing temperature. We also show the charge density N(t) 
in Fig. 7: solid line (ksT = 25 meV) and dashed line 
(k B T = 10 meV). We see that N(t) = N L (t) + N R (t) 
approaches the same value for fcsT = 10 meV and 
ksT = 25 meV in the steady state. To understand this 
feature, we plot Nl and Nr, which are defined as the 
charge density provided from the left and right electrodes, 
respectively. In the steady state (or a large time) , Nl de- 
clines with increasing temperature, on the other hand Nr 
increases with increasing temperature. We see that they 
compensate in each other. 

Next, we study the transient current for the case of 
different tunneling rate ratio T l/ Gamma r at room tem- 
perature ksT = 25 meV. Comparing with the results at 
zero temperature, the oscillation structure of J(t) at room 
temperature is more complicate due to the enhancement 
of Jj?(i). For T L = 0.7 meV and T L = 0.3 meV , the for- 
mer exhibits the exponential decay, the latter exhibits the 
exponential growth. This indicates that the location of 
QD exists a considerable effect on the transient transport 
properties of SETs. Even though the above theoretical 
analysis is for SETs, we attempt to report the experimen- 
tal measurement of the transient current of Ge SHT. For 
multivalley conduction band of germanium semiconduc- 
tors, the intervalley interaction effect can not be ignored. 
However, such effect is not included in Eq. (1). The 
more realistic system described by Eq. (1) is a Ge SHT, 
where valence band is a single valley. In Figure 9 we show 
the transient current of Ge SHT for two different applied 
voltages V a = 20 mV and V a = 30 mV at room temper- 
ature and V g — — 3 V. The detailed fabrication process 
of Ge SHTs and measurement technique of transient cur- 
rent will be discussed in elsewhere 7 . The solid lines are 
the experimental curves. The dashed lines are the fitting 
curves. According to fitting curves, the transient currents 
display exponential growth feature. Based on previous 
theoretical analysis, the Ge SHT has the characteristic of 
Fl/Fr < 1- Two fitting curves also provide the bias- 
dependent tunneling times 70 sec and 60 sec for the ap- 
plied bias V a = 20 mV and V a = 30 mV. This very long 
tunneling time is arising from a fact of very high barrier 
of SiO-2 which the Ge QD is embedded into. Besides, we 
observe the small oscillation structures on the exponen- 



tial growth curves. The oscillation structures vanish for 
sufficient long time. This indicates that this oscillation 
structures are not the measurement error. Owing to the 
lack of information about the detailed size and shape of 
QD, there are some difficulties to do the detailed compar- 
ison between the theoretical calculation and experimental 
measurement. 

IV. SUMMARY 

In this study the transient current of a SET is investi- 
gated by using the formula derived by Jauho, Wingrcen 
and Meir. 2 The tunneling rate ratio exists a considerable 
influence on the transient current of SETs. The tunneling 
current of a single Ge SHT at room temperature displays 
the oscillation structures in the transient process. This 
feature is attributed to time-dependent density of states 
of QD in a switching on the gate voltage. 
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Figure Captions 

Fig. 1. Transient current through the symmetry bar- 
rier for different applied voltages at zero temperature and 
V g = 25 mV: solid line (V a = 30 mV), dashed line 
(V a = 25 mV) and dotted line (V a = 20 mV), where 
J = e meV/h and to = V (meV). 

Fig. 2. Transient current through the symmetry barrier 
for V g = 25 mV, V a — 30 mV and zero temperature. 
Solid line is duplicated from that of Fig. 1. Dashed line 
and dotted line denote J L \(t) and J L \. Dash-dotted line 
represents Jjp(i). 
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Fig. 3. Transient current through the symmetry barrier 
for different applied gate voltages at zero temperature and 
V a = 30 mV: solid line (V g = 25 mV), dashed line (V g = 
20 mV), dotted line (V g — 15 mV) and dash-dotted line 
(V g = 10 mV), where Jo/e = meV/Ti and to = h/(meV). 

Fig. 4. Transient electron occupation number for 
different applied gate voltages at zero temperature and 
V a = 30 mV. The curves are one to one corresponding to 
those curves of Fig. 2. 

Fig. 5. Transient current for different tunneling rate 
ratio at zero temperature, V a — 30 mV and V g = 25'mV: 
the curves from the bottom to the top (for t/to < 1) cor- 
respond, respectively, Tl = 0.3,0.4,0.5,0.6 and 0.7 meV. 
Meanwhile, r = T L +T R = 1 meV. Here J /e = meV/h 
and to = h/(meV). 

Fig. 6. Transient current through the symmetry 
barrier for different temperatures at V a — 30 mV and 
V g = 25 mV: solid line (fc^T = meV), dashed line 
(k B T = 15 meV), dotted line (k B T = 20 meV) and dash- 
dotted line (k B T = 25 meV), where Jo/e — meV/h and 
t = h/(meV). 

Fig. 7. Transient electron occupation number for two 
different temperatures: solid line (k B T — 25 meV) and 
dashed line (k B T = 10 meV). 

Fig. 8. Transient current for different tunneling rate 
ratio at temperature k B T — 25 meV, V a = 30 mV 
and V g = 25'mV. T = T L +T R = 1 meV . Here 
Jo/e = meV/h and t — h/(meV). 

Fig. 9. Transient current of Ge single hole transistor for 
different applied voltages at room temperature and V g = 
3 V. The curves from the bottom to the top correspond 
to V a = 20 mV and 30 mV, respectively. 
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